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Accurate measurement of the frequency -dependent shift of the self-absorbed radio core is required 
for multi-frequency analysis of VLBI data since absolute positional information is lost as a result 
of phase self-calibration. We use the cross-correlation technique of Croke & Gabuzda (2008) 
on the optically thin jet emission to align our VLBA images. Our results are consistent with 
those obtained from the phase-referencing method, as well as alignment by model-fitted optically 
thin jet components. Physical parameters of the compact jet regions, such as the magnetic field 
strength (B) and the distance of the radio core to the jet origin (r), can be calculated from these 
measurements. For the source Mrk 501, we find a magnetic field strength of 0.15 ±0.04 G in 
the 8.4-GHz core at a distance of 0.8 ± 0.2 pc from the base of the jet. By extrapolating our 4.6 
to 15.4 GHz results for BL Lac (2200+420), we estimate magnetic field strengths of the order of 
1 G in the millimetre VLBI core. Using our core-shift measurement between 1.6 and 4.8 GHz for 
1803+784, we find fi C ore(4.8 GHz) = 0.11 ±0.02 G and r eore (4.8 GHz) = 20±5 pc. The phase- 
referencing observations of this source at 8.4 and 43 GHz by Jimenez -Monferrer et al. (2008) 
imply B cor e(43 GHz) = 1.0 ±0.4 G and r core (43 GHz) = 2.0 ±0.9 pc. 
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1. Introduction 

Highly collimated jets of relativistic plasma are considered to be generated and propelled out- 
wards from the central regions of Active Galactic Nuclei (AGN) by magnetic forces surrounding 
the supermassive black hole, see [1] and references therein. High resolution VLBI observations 
can detect the synchrotron emission from these jets in the radio regime at distances of the or- 
der of 10 3 — 10 7 Rs, where R$ is the Schwarzchild radius of the black hole, corresponding to the 
sub-parsec to parsec scales of the jet [2]. Current theoretical models [3, 4] invoke rotation of the 
magnetosphere around the compact central object which creates a stiff helical magnetic field that 
expels and collimates the jet flow. At a distance of several hundred gravitational radii from the cen- 
tral engine, the natural pinching action of the magnetic field can over-collimate the flow resulting 
in a "collimation shock" which can either partially or completely disrupt the helical field structure 
in the jet. Recent observational evidence suggests that this "collimation shock" may occur in the 
millimetre-wave VLBI core [5]. A significant number of VLBI observations at cm-wavelenghts 
support the presence of helical magnetic fields in the parsec scale radio jet (e.g., [6] and refer- 
ences therein), indicating that the initially ordered magnetic field is not completely disrupted and 
that either remnants of the earlier field structure remain or possibly a current driven helical kink 
instability is generated [7-9]. 

The radio cores of AGN are generally flat-spectrum, self-absorbed regions. Since VLBI res- 
olution is usually not sufficient to completely resolve the true optically thick radio core, the core 
region contains emission from around the optical depth (t) = 1 surface and some contribution 
from the optically thin inner jet. From the jet model of [10], the frequency dependence of the 
self-absorbed core is described by r oc y~ l / k ' ', where r is the distance from the central engine and 
k r = ((3 — 2a)m + 2n — 2) / (5 — 2a) with m and n describing the power-law fall-off in the magnetic 
field strength and particle number density, respectively, with distance from the central engine. For 
equipartition between the jet particle and magnetic field energy densities, the quantity k r = l, with 
the choice of m = 1 and n = 2 being reasonable (e.g., [11 - 13]), making k r independent of the 
spectral index. Therefore, in equipartition, the core-shift (Ar — \f~vi — ^V?|) between two frequen- 
cies (V2 > Vi) is directly proportional to (V2 — Vi)/(V2Vi). Hence, for accurate analysis of multi- 
frequency VLBI maps, we must correct for the frequency-dependent position of the self-absorbed 
radio core. 

2. Alignment Method & Results 

The frequency dependent core-shift is generally obtained by aligning model-fitted optically 
thin jet components [13], but this can be rather difficult for jets that lack distinct components across 
all frequencies. Phase -reference observations [14] can be used to obtain the absolute position of 
the radio core by comparing the position of the target source with respect to a reference source, 
which should ideally be a point source, but this is a very calibration-intensive approach. 

Aligning the images via the positions of optically thin jet components can be very difficult due 
to the complicated total intensity structure of VLBI jets, especially at low frequencies where often 
there are no obvious, common optically thin features in the diffuse emission region. However, 
using a cross-correlation technique for aligning the optically thin emission from the jet, see [15] 
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for details, we have been able to find reliable measurements for the frequency dependent core-shift 
of a number of sources across a wide range of frequencies. 

2.1 Mrk501 

We have found core-shift measurements for the well studied source Mrk 501 (1652+398) from 
multi-frequency VLBA observations [16] obtained in May 1998. For this source, we are able to 
compare our alignment results with those based on the alignment of individual model-fitted jet 
components. Using space- VLB I observations of Mrk 501 at 1.6 GHz and ground-based observa- 
tions at 5 GHz in April 1998, [17] were able to align these two images by aligning the positions 
of the optically thin jet components. This was possible due to the similarity in resolution of the 
5 and 1.6 GHz images provided by the space baselines at 1.6 GHz. Our observations of Mrk 501 
at 1 .6 and 5 GHz could not be aligned in this manner due to the lack of distinct jet components 
in the poorer resolution 1.6 GHz image. Instead, we used the cross-correlation technique to align 
the images and found a core-shift of 0.78 ± 0.07 mas. This value compares very well with the 
measured core-shift of 0.72 ± 0.28 mas from the aligned images of [17]. 

2.2 BL Lac 

BL Lac (2200+420) has a smooth, extended jet structure lacking a bright distinct jet knot, 
across a wide range of frequencies (e.g., [6]), with which to align the images using the model-fitted 
component method. However, the VLBA jet images of this source are straightforwardly aligned 
using the cross-correlation technique. We have analysed the core-shift using VLBA observations 
at 4.6, 5.1, 7.9, 8.9, 12.9 & 15.4 GHz. The excellent frequency coverage provides a degeneracy 
of core-shift measurements with different frequency combinations reinforcing each other. The 
magnitude of the separation between the 4.6 and 15.4 GHz core is 0.44 ±0.07 mas. The direction 
of the shift is 16.3 ± 0.9° NE, which is consistent with the inner jet direction at 43 GHz [6] of 
197 ± 1° (180° + 17°) at this epoch. 

2.3 1803+784 

Using VLBA observations of the AGN 1803+784 at 1.6 and 5 GHz, we derived a core-shift 
of 1.27 ±0.35 mas using the cross-correlation technique. We are able to indirectly test the ac- 
curacy of our result from independent phase-referencing observations for this source [18]. They 
obtained a core-shift of 0.27 ±0.13 mas between 8.4 and 43 GHz from their astrometric analysis 
for this source. Using our lower frequency measurement and assuming equipartition between the 
jet particle and magnetic field energy densities, we predict a shift of 0.29 mas between 8.4 and 43 
GHz, consistent with the phase-referencing results. Conversely, this provides some support for the 
assumption that equipartition holds for the compact jet region of 1803+784 between 1.6 and 43 
GHz. However, it should be noted that the value of the core-shift may change at different epochs 
due to nuclear flares [19]. 

3. Jet Physics 

We use the equations described in [13], from the model of [11], in order to calculate some 
physical parameters of the jets listed above using the core-shift measurement. Knowledge of the 
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Figure 1: Plot of core-shift versus frequency 
for Mrk 501 using 8.4 GHz as the reference fre- 
quency. Red line: Ar — A(v~ l l kr — 8.4 _1 / A> ) with 
best-fit parameters of A = 1.98 ±0.21 and k r = 
0.82±0.11 



Figure 2: Plot of core-shift versus frequency 
for BL Lac using 15.4 GHz as the reference fre- 
quency. Red line: Ar = A(v~ 1//A> — \5A~ l l kr ) 
with best-fit parameters of A = 2.97 ±0.19 and 
k r = 0.98 ±0.06 



redshift of the source as well as estimates of the jet opening angle, viewing angle and bulk Lorentz 
factor are required to obtain estimates of the magnetic field strength and other quantities such as 
the distance of the radio core to the base of the jet. 

The parameter k r , described in Section 1, can be calculated from core-shift measurements at 
several different frequencies. Therefore, instead of assuming the equipartition value of k r = 1, we 
can measure the exact value. In the case of Mrk 501, we analysed the core-shift from observations 
at 1.6, 2.2, 4.8 and 8.4 GHz. Using 8.4 GHz as the reference frequency, we obtain a value of 
k r = 0.82 ±0.11 (Figure 1). This result is consistent with the synchrotron self-compton (SSC) 
models of [20], who find best fit values of m = 0.9 and n = 1.8 for the jet of Mrk 501, which 
gives a value of k r = 0.86 using an optically thin jet spectral index of a = —0.7. Using the well 
constrained jet parameters from [17], we find a magnetic field strength of 0.15 ± 0.04 G in the 
8.4-GHz core at a de-projected distance of 0.8 ± 0.2 pc from the jet origin. 

We perform a similar analysis for BL Lac using observations from 4.6 to 15.4 GHz. In this 
case, we find a value of k r = 0.98 ± 0.06 consistent with the equipartition value of k r = 1 (Figure 
2). Using the jet parameters from [21], we estimate the magnetic field strength in the VLBI core 
at each frequency and its distance from the base of the jet. Figure 3 shows a plot of the magnetic 
field strength at various distances along the jet. Extrapolating this plot, we find that the magnetic 
field strength approaches ~1 G in the mm-core. This is consistent with results from [22], where he 
finds magnetic field strengths of order 1 G in the mm-core of 3C 273 by measuring the synchrotron 
self- absorption (SSA) turnover frequency of individual model-fitted components [23]. Importantly, 
this method does not require the assumption of equipartition. 

Even though we have not found a value of k r for 1803+784, our results in Section 2.3 suggest 
that an equipartition value of k r = 1 is not unreasonable. Hence, we find B core (4.8 GHz) = 0.11 ± 
0.02 G and r core (4.8 GHz) = 20 ±5 pc using jet parameters from [24]. The phase-referencing 
measurement of [18] implies a magnetic field strength in the 43-GHz core of 1.0 ±0.4 G at a 
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15.4 GHz 8.9 GHz 5.1 GHz 
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Figure 3: Plot of magnetic field strength (B) in units 
of Gauss versus distance along the jet (r) in units of 
parsecs for the core of BL Lac at 5.1, 7.9, 8.9, 12.9 
and 15.4 GHz. Also plotted is the line B = B\ ?Q r~ x 
with B lpc = 0.4 G. 
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Figure 4: Straight-line fit to the core-shift 
data for Mrk 501 assuming equipartition. The y- 
intercept value is equal to —0.03 ± 0.02 mas. A 
straight-line fit through the origin is expected for 
k r exactly equal to 1 . 



distance of 2.0 ± 0.9 pc from the base of the jet. 

One of the limitations of the core-shift method is the need to assume equipartition in order to 
find an estimate for the particle number density and hence, the magnetic field strength. In the case 
of Mrk 501, the value of k r doesn't deviate too strongly from the equipartition value (see Figure 
4) but, for this source, estimating the magnetic field strength using the SSA turnover method may 
provide more reliable estimates. At least in the case of 2200+420, we know that the compact jet 
region is in the equipartition regime between 4.6 and 15.4 GHz; therefore, we have confidence 
that our magnetic field estimates are accurate. For 1803+784, by combining our measurements 
with phase-referencing observations, we have indirect evidence that the compact jet region does 
not deviate strongly from the equipartition regime between 1.6 and 43 GHz. 



4. Conclusions 



Using a cross-correlation technique for aligning the optically-thin jet emission, we obtain ac- 
curate measurements of the frequency-dependent shift of the self-absorbed VLBI core of three 
AGN jets. We have tested our method against both phase-referencing observations and alignment 
by individual model-fitted steep-spectrum jet components. Using our measurements of the core- 
shift, we find magnetic field strengths in the compact centimetre jet regions of Mrk 501, BL Lac 
and 1803+784 of the order of 100's of mG; with values approaching 1 G in the mm-wave VLBI 
cores of BL Lac and 1803+784. We find that equipartition is a valid assumption for the compact 
inner jet regions of BL Lac, while we find some deviation from equipartition in the case of Mrk 
501. We have also found estimates for distances of the observed VLBI core to the base of the ra- 
dio jet for each source, providing the approximate location of the central supermassive black hole. 
Estimates such as these are essential in efforts to test the predictions of current theoretical models 
on the locations of regions of jet launching, acceleration and collimation. 
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